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Nitrogen- Vacancy color centers in diamond 
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We report a study of the ''E excited-state structure of single negatively-charged nitrogen- vacancy 
(NV) defects in diamond, combining resonant excitation at cryogenic temperatures and optically 
detected magnetic resonance. A theoretical model is developed and shows excellent agreement with 
experimental observations. Besides, we show that the two orbital branches associated with the '^E 
excited-state are averaged when operating at room temperature. This study leads to an improved 
physical understanding of the NV defect electronic structure, which is invaluable for the development 
of diamond-based quantum information processing. 
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Coupling between flying and stationary qubits is one 
of the crucial requirements for scalable quantum infor- 
mation processing [Jl |5] . Among many quantum systems 
including single atoms [3^ and semiconductor quantum 
dots [4J, the negatively-charged nitrogen- vacancy (NV) 
color center in diamond is a promising solid-state candi- 
date for realizing such interface, owing to long spin co- 
herence time of their spin states [5j and availability of a 
strong optical transition [6J. Nevertheless, even though 
NV defects have been intensively studied during the last 
decades, the excited-state structure as well as the dynam- 
ics of excitation-emission cycles are surprisingly not yet 
fully understood. This knowledge is however of crucial 
importance for the realization of long-distance entangle- 
ment protocols based on coupling of spin-state to optical 
transitions [HEIE]- 

In this Letter, we report a study of the excited-state 
structure of single NV defects as a function of local strain, 
combining resonant excitation at cryogenic temperatures 
and optically detected magnetic resonance (ODMR). Be- 
sides, we show that the two orbital branches associated 
with the excited-state are averaged at room tempera- 
ture. A theoretical model is developed and shows excel- 
lent correspondence with experimental observations. 

The NV color center in diamond consists of a substi- 
tutional nitrogen atom (N) associated with a vacancy 
(V) in an adjacent lattice site, giving a defect with Csv 
symmetry. For the negatively-charged NV color center 
addressed in this study, the ground state is a spin triplet 
■^A2 [inj [HJ [12] . Spin-spin interaction splits ground state 
spin sublevels by 2.88 GHz into a spin singlet Sz, where 
z corresponds to the NV symmetry axis, and a spin dou- 
blet Sx,Sy (see Fig. [ija)). The excited state '^E is also 
a spin triplet, associated with a broadband photolumi- 
nescence emission with zero phonon line (ZPL) around 
637 nm (1.945 eV). Besides, the '^E excited state is an 
orbital doublet, which degeneracy is lifted by non-axial 
strain into two orbital branches, and Ey, each or- 



bital branch being formed by three spin states Sx, Sy and 
Sz (see Fig. [l|a)) [13]. As optical transitions ^A2 ^^E 
are spin-conserving, excitation spectra of single NV color 
centers might show six resonant lines, corresponding to 
transitions between identical spin sublevels. 

The order of other energy levels is still under debate 
but it is now well established that at least one metastable 
state ^Ai is lying between the ground and excited triplet 
states [11]. Non-radiative inter-system crossing to the 
^Ai state is strongly spin selective as the shelving rate 
from the Sz sublevel is much smaller than those from Sx 
and Sy. Furthermore, the metastable state decays pref- 
erentially towards the ground state spin level Sz , leading 
to a strong spin polarization into that state after a few 
optical excitation-emission cycles. Thereby, optical tran- 
sitions linking Sx or Sy sublevels are non cycling tran- 
sitions. Experimental investigation of the excited-state 
structure then requires the use of a microwave excitation 
resonant with the ground state transition at 2.88 GHz, 
in order to maintain a non-zero time-averaged population 
within each of the ground state spin sublevels [13J. 

Spectroscopic investigations require good spectral sta- 
bility of NV defects. Using low nitrogen concentra- 
tion type Ila natural diamond, Fourier-transform lim- 
ited emission have been recently reported at cryogenic 
temperature [151 I16j . In the present study, we investi- 
gate native single NV defects in an ultra-pure synthetic 
type Ila diamond crystal prepared using a microwave as- 
sisted chemical vapor deposition process. Within such 
sample, where the nitrogen concentration is below Ippb 
(< 1.7 X lO^^cm"'^) and the NV defect concentration 
smaller than 10^°cm~^, all single NV defects have shown 
perfect spectral stability. 

NV defects are addressed using confocal microscopy at 
cryogenic temperatures (T w 4 K) . A tunable laser diode 
is used to excite resonantly the NV centers on their ZPL. 
The red-shifted photoluminescence (PL) between 650 nm 
and 750 nm is detected in a confocal arrangement and 



2 




0.0 0.1 0.2 0.0 0.1 0.2 

MW pulse duration (us) MW pulse duration (|is) 



FIG. 1; (color online). (a)-Energy- level diagram of a single 
NV defect. (b)-Excitation spectrum of a single NV defect 
presenting six resonances, which correspond to optical transi- 
tions linking identical spin sublevels. Solid line is data fitting 
using Lorentzian functions. Energy splittings between each 
excited-state spin sublevels for this NV center are written in 
blue in (a) (not in scale). (c)-Pulse sequence used to measure 
single-electron spin Rabi oscillations (see main text), (d) and 
(e)-Rabi nutations measured using resonant optical transi- 
tions for the spin read-out in the upper branch (d) and in 
the lower branch Ey (e). A tt phase shift of the Rabi nutation 
allows to discriminate for each orbital branch between optical 
transitions linking Sz or Sx,y spin sublevels. 



used to monitor excitation spectra by sweeping the laser 
diode frequency. In addition, microwaves (MW) reso- 
nant with the ground state spin transition are appHed 
via a copper microwire located close to the NV defects. 

A typical excitation spectrum of a single NV center is 
depicted in Fig. [TJb). As expected, three transitions are 
observed for each excited-state orbital and Ey, corre- 
sponding to transitions between identical spin sublevels. 
Unambiguous assignment of optical transitions linking 
Sz or Sx,Sy spin sublevels is made by measuring single 
electron spin Rabi oscillations using the pulse sequence 
shown in Fig. [l|c). The NV center is first initialized into 
the ground state Sz sublevel using an optical pulse of du- 
ration 3 fis at the wavelength A = 532 nm. A microwave 
pulse at 2.88 GHz is then applied and the spin-state is 
finally read-out by measuring the PL intensity using a 
1 /iS red laser pulse resonant with a given transition. De- 
pending on the optical transion used for the spin read-out 
{Sz or Sx, Sy), the Rabi nutation shows a tt phase shift, 
allowing to discriminate between optical transitions link- 
ing different spin sublevels (see Fig. [ijd)-(e)). 



The energy level scheme presented in Fig. [TJa) is only 
valid for a particular local strain of the diamond lattice 
in the vicinity of the NV defect. Indeed, non-axial local 
strain strongly modifies the energy difference between the 
two orbitals E,j; and Ey , as well as the relative position of 
the spin sublevels. Following a model previously devel- 
oped in Ref . [T3] , we investigate the position of the energy 
levels as a function of the local strain. The Hamiltonian 
of the '^E excited-state is given by 



H — Hq + Hso + Hss + Hstr 



(1) 



where Hq is the dominant term giving the energy of 1.945 
eV with respect to the ''A2 ground state, Hso is the spin- 
orbit coupling, Hss the spin-spin interaction, and Hstr 
the perturbation resulting from local strain. The axial 
spin-orbit interaction XzLzSz sphts the excited-state 
into three twofold degenerate levels, noted {£), {£') and 
{Ai,A2) according to the irreducible representations of 
Csy symmetry (see Fig. [2]). The transverse spin-orbit 
interaction X^.y is weak {Xx.y = 0.2 GHz |13j ) and can 
be neglected at present. Spin-spin interactions displace 
all the latter levels in a similar fashion as in the ground 
state by the energy D^g{Sl — 2/3). Besides H^^ lifts the 
degeneracy of the Ai and A2 states by ± A [TS] . The fine 
structure at zero strain is therefore given by the three 
parameters Xz, D^s and A (see Fig. [2]). 

The effect of local strain can be described using the 
Hamiltonian Hgtr = S^V^ + SyVy + SzVz, where Si and 
Vi represent respectively the strain parameter and an or- 
bital operator in the i direction [18]. Axial strain Sz gives 
rise to a linear shift of all energy levels and is not con- 
sidered further. Transverse strains, Sx and Sy, split the 
energy levels into two spin triplets, associated with the 
two orbital branches Ex and Ey. The energy splitting is 
given by ±i5j_, where S± cx {S^ + Sy)i (see Fig. ^ [T5]. 

Excitation spectra of several single NV defects have 
been measured, corresponding to different values of the 
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FIG. 2: (color online). Schematic diagram of the excited-state 
energy levels including spin-orbit interaction, spin-spin inter- 
action and transverse local strain. Note that 1 GPa external 
stress gives approximately 10"^ GHz splitting [TT]. 
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FIG. 3; (color online). (a)-Excitation spectra recorded for 
different single NV defects denoted NVl, NV2 and NV3. (b)- 
Excited-state energy structure as a function of the transverse 
strain S± for 27 single NV defects. Solid lines correspond to 
the calculation without any free parameter and the inset is 
a zoom on the lower branch for a strain corresponding to an 
avoided crossing between \Ey, Sz) and \Ey,Sa:} spin sublevels. 
At this avoided crossing, the solid lines exchange colors for 
clarity purpose. (c)-Energy level scheme of the lower branch 
for the defect NV3, showing three spin conserving transitions 
(solid arrows) and two spin- flip transitions (dashed arrows), 
appearing when the states \Ey, Sz) and \Ey, Sx) are crossing. 
(d)-Energy level scheme of the upper branch for the defect 
NV2. The energy splitting between the states \Ex,Sz) and 
\Ej:,Sy) is equal to 2.88 GHz, leading to a repumping effect 
of the \Sz) — > \Ex,Sz) optical transition. 



transverse strain S± (see Fig. [3]^a)-(b)). From low strain 
experimental data, we estimate the zero-strain parame- 
ters to be = 5.3 GHz, D^^ = 1-42 GHz and A = 
1.55 GHz. Using such measurements, a calculation of 
the eigenenergies of the Hamiltonian described by equa- 
tion ([T| is performed for different strengths of the trans- 
verse strain 6±, including the effect of transverse spin- 
orbit interaction X^ y = 0.2 GHz |T3]. All measurements 
are in excellent agreement with this calculation, as de- 
picted in figure [3]jb) . 

With such results, we have a complete physical under- 
standing of the NV defect dynamics in terms of cycling 
spin conserving transitions and spin- flip transitions. 

In the lower branch Ey, spin sublevels are strongly 
mixed by local strain, resulting in non-cycling transi- 
tions, as spin-flip can occur either by radiative decay 
or by non-radiative decay through the metastable state. 
Thus, the lower branch cannot be observed in excita- 
tion spectra without applying MW. Moreover, for a par- 
ticular strain corresponding to an avoided crossing be- 
tween \Ey,Sz) and \Ey,Sx) spin sublevels (see inset of 



Fig.|3jb)), the mixing becomes so important that fully al- 
lowed spin-flip optical transitions can be observed. This 
situation corresponds to the excitation spectrum mea- 
sured for the defect NV3 (see Fig.|3][a) and (c)) where two 
additional spin-flip transitions are observed in the lower 
branch. We note that the same effect occurs at higher 
strain when \Ey, Sz) and \Ey, Sy) spin sublevels are cross- 
ing. Such spin-flip transitions can be used for single- 
spin high-speed coherent optical manipulation through 
A-based scheme HHHD]. 

In the upper branch E^, the situation is completely 
different as no crossing occurs between different spin- 
sublevels. When the MW is not applied, a single res- 
onant line remains visible in the excitation spectrum, 
corresponding to the cycling transition \Sz) \Ex,Sz). 
Nonetheless, for most of the studied NV defects this reso- 
nance is observed to be several times weaker without ap- 
plying MW. This is explained by a remaining small mix- 
ing of the excited-state spin sublevels through non-axial 
spin orbit interaction X^^y |10j . Remarkably, for a spe- 
cific local strain, the optical transition l^^) \Ex,Sz) 
is however found to be much stronger without applying 
MW. This situation is described in Fig.|3|d), correspond- 
ing to the excitation spectrum of defect NV2 shown in 
Fig. [3]ja). At this particular strain the energy splitting 
between \Ex,Sz) and \Ex,Sy) is equal to 2.88 GHz, like 
in the ground state. In such a configuration, if a spin-flip 
occurs during \Sz) — > \Ex,Sz) transitions, the NV defect 
is efficiently repumped in the Sz ground state through a 
resonant \Sy) \Ex,Sy) optical transition followed by 
non-radiative decay to the metastable state responsible 
for spin-polarization of the NV defect. Note that local 
strain can be tuned by applying electric field through 
electrodes deposited on the diamond surface, allowing to 
externally control conditions of spin-conserving cycling 
transitions or spin- flip transitions for A scheme |13j . 

Excited-state spectroscopy of single NV defect using 
ODMR techniques at room temperature has recently un- 
derlined an excited-state electron spin resonance (ESR) 
around 1.4 GHz f2T1'22|. Surprisingly, low temperature 
excitation spectra described above do not show any sig- 
nature of such an energy splitting between excited-state 
spin sublevels, neither in the upper branch nor in the 
lower branch. In order to understand this observation, 
ODMR spectra are measured for the same single NV 
defect at different temperatures by sweeping the MW 
frequency. As depicted in Fig. Qa), when the experi- 
ment is performed close to room temperature, the well- 
known ground state ESR is detected at 2.88 GHz, and 
the excited-state ESR around 1.4 GHz is observed, as 
previously reported. Performing the same experiment at 
low temperature (T= 6 K) shows that the excited-state 
ESR fully disappears. The evolution of the excited-state 
ESR contrast as a function of temperature is depicted in 
Fig.Qb), showing a strong decrease around T = 150 K. 

This observations are explained by considering that the 
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FIG. 4: (color online). (a)-ODMR spectra of a single NV color 
center recorded at T = 260 K and at T = 6 K. (b)-Contrast 
of the excited-state ESR as a function of temperature. The 
contrast is normalized to the integral of the ground state ESR 
at 2.88 GHz. (c)-Energy splitting between the averaged en- 
ergy of Sz and Sx,y sublevels of the two orbital branches as 
a function of strain. Data points are inferred using the mea- 
surements depicted in Fig.jsjb). The black dashed line corre- 
sponds to Des = 1.42 GIfz. 



The physical process responsible for orbital averaging, 
which remains under question, will be addressed in future 
work. This process, which could induce flips between the 
two orbitals at a frequency much higher than the radia- 
tive lifetime without altering the spin projection, could 
be for example a phonon-mediated process [53]. Similar 
process would also account for why ESR is not observed 
in the gound state of the neutral NV defect [24] . 

Summarizing, we have reported a study of the 
negatively-charged NV defect excited-state structure as 
a function of local strain. The present work gives sig- 
nificant insights into the NV defect electronic structure, 
which is invaluable for the development of diamond-based 
quantum information processing. 
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(SFB/TR21) and Australian Research Council. V. J. ac- 
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two excited-state branches are averaged when operating 
at room temperature [23] . Room temperature studies 
have indeed shown that the isotropic g-factor associated 
with the excited-state ESR at 1.4 GHz is similar to the 
ground state g-factor {g « 2) [^TJ [22], indicating that 
the orbital angular momentum does not play a signifi- 
cant role in the excited state at room temperature. We 
now assess the expected value of the energy splitting be- 
tween Sz and Sx , Sy excited-state spin sublevels by av- 
eraging over the orbital branches. Such an orbital av- 
eraging quenches both the effect of spin-orbit splitting 
A2 and the effect of spin-spin interaction A, which splits 
Ai and A2 sublevels [TH] (see Fig. [2|. The averaged en- 
ergy splitting between Sz and Sx,Sy excited-state spin 
sublevels is then only determined by the parameter Des, 
which has a value of 1.42 GHz. This averaged value can 
be checked using the measurements shown in Fig. [sjb), 
by estimating the energy splitting between averaged Sz 
and Sx, Sy spin sublevels. Such splitting is found around 
1.4 GHz for the whole range of considered local strains, 
as depicted in Fig.|4|^c). This is consistent with the value 
of Des and, hence, the averaging process explains why an 
excited-state ESR is detected around 1.4 GHz at room 
temperature. We note that for very large strain, which 
leads to a lowering of Ca^ symmetry, the spin-spin inter- 
action becomes Hss = Des{S'^—2/3)+Ees{Sx — Sy) where 
Ees is proportional to the transverse strain. In such con- 
ditions, the averaged energy splitting is then given by 
Des ± -Eg. In case of large strain, the orbital averaging 
then accounts for previously reported observations of a 
splitted excited-state ESR at room temperature [HJ [52] . 
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